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In a recent manuscript (arXiv: 1208.5046) Peter Sorensen claims that XENONlOO's upper limits 
on spin-independent WIMP-micleon cross sections for WIMP masses below 10 GeV "may be un- 
derstated by one order of magnitude or more" . Having performed a similar, though more detailed 
analysis prior to the submission of our new result (arXiv: 1207.5988), we do not confirm these find- 
ings. We point out the rationale for not considering the described effect in our final analysis and 
list several potential problems with his study. 
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In a recent manuscript [1], P. Sorensen examines 
our results from a 225 live-days dark matter run with 
XENON100 [2] and claims that the XENON100 upper 
limit on WIMP-nucleon cross sections at WIMP masses 
below 10 GeV might be significantly stronger than our 
published result. We are aware of the raised issues and 
take the opportunity to comment here. While we wel- 
come the author's endorsement of our main conclusion, 
namely the lack of an observed dark matter signal in this 
run, we do not support his statement of one order of mag- 
nitude improvement in sensitivity for low-mass WIMPs 
after having performed a similar analysis prior to the 
submission of our manuscript to PRL. 

We agree with the argument that in principle one 
might use the additional information carried by the pro- 
portional light signal, S2, in order to obtain a better mea- 
sure of the energy of each scattering event in our detector. 
We would thus exploit not only SI, the prompt scintilla- 



tion signal, but the fully available phase space. Indeed, 
as shown in [3] we have used the combined SI and S2 in- 
formation to significantly improve the energy resolution 
of our detector for interactions of gamma rays at various 
energies and to understand its main background sources 
[4]. On the other hand, as we discuss in more detail 
later, we are still unable to use the information in the 
S2-channel at the energies of interest to a dark matter 
search, for we lack measurements of the ionization yield, 
Q y , of liquid xenon for nuclear recoils of a few keV. We 
also agree with the statement that low-mass WIMPs are 
expected to show a different S2/S1 versus SI distribu- 
tion than the one expected from calibration data with an 
241 AmBe neutron source. In fact, we have studied these 
effects in detail, in a similar fashion as followed in the 
paper by P. Sorensen: we have inferred Q y based on our 
241 AmBe nuclear recoil calibration data and on the mea- 
sured C c f{(E nr: ) and have used Monte Carlo simulations to 
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generate the expected distributions of events in the S2/S1 
versus SI parameter space for various WIMP masses. 
We have analyzed the influence on our upper limit on 
WIMP-nucleon cross sections using the same profile like- 
lihood method as described in [2, 5]. For a WIMP mass of 
6GeV, the upper limit changes from 8.33x 10 _40 cm 2 to 
6.45x 10 _40 cm 2 , at 90% confidence level, which means an 
improvement by a factor of 1.3. The difference becomes 
smaller with increasing WIMP masses, being a factor of 
1.2 for a 50GeV WIMP. 

At present, we do not consider this effect in our final 
analysis reported in [2] for the following reasons: 

• No direct measurements of Q y at the low nuclear 
recoil energies relevant for the low-mass WIMP re- 
gion exist, as also shown in Figure 2 (upper panel) 
in the P. Sorensen manuscript. We thus depend on 
its inference based on C c g and the measured de- 
tector response to an 241 AmBe neutron calibration 
source. Such an inference, while a priori viable, 
requires the simultaneous disentanglement of the 
Q y behaviour at low energies from efficiencies in 
SI and S2 and thus introduces additional system- 
atic errors into the analysis. We remark that direct 
measurements of Q yi using small liquid xenon time 
projection chambers (TPCs) are in progress within 
our collaboration. Such an independent measure- 
ment of the ionization yield will greatly reduce the 
systematics involved in this approach and will make 
it a viable option for future analyses of dark matter 
data from liquid xenon TPCs. In contrast to Q y , 
direct measurements of the light yield, parameter- 
ized via £ c g, exist [6], the lowest measured data 
point being at 3 keV nuclear recoil energy [7] . In- 
direct determinations of £ e ff, as attempted in [8] 
and [9], have not yielded results compatible with 
direct measurements [7]. These difficulties are fur- 
ther high-lighted by a re-analysis of the result in [9] 
for the same experiment [10], which are now com- 
patible with data from direct C e g determinations. 

• The reason why the 241 AmBe calibration data do 
not accurately describe the expected S2/S1 versus 
SI distribution of low-mass WIMPs can be mainly 
ascribed to four factors: i) The threshold in SI 
is much higher in energy than the corresponding 
threshold in S2; ii) The resolution in SI in our 
TPC is much coarser than the resolution in S2; iii) 
The steep energy spectra and kinematic cut-off of 
low-mass WIMPs are well below the mean energy 
threshold set by the SI measurement; iv) S2 signals 
do not fluctuate upwards coherently with SI, but 
are rather independent processes. As a result, an 
upward fluctuation of a nuclear recoil well below 
the mean energy threshold in SI will most likely 
result in an S2 close to its mean corresponding to 
the low energy. The corresponding threshold in S2 
is much lower, allowing us to sample the bulk of 
the distribution in S2 but only the tail in SI. The 



response for low-mass WIMPs therefore systemat- 
ically moves towards the lower left corner of our 
two-dimensional dataspace. 

• We have performed a profile likelihood analysis [5] 
for the case of WIMP-induced interactions dis- 
tributed qualitatively similar to Figure 4 in [1]. 
In this study, we have used a mean C e g(E nl ) as 
shown in [2] and a Q y derived from comparing our 
241 AmBe calibration data to a dedicated Monte 
Carlo study which includes detailed knowledge of 
our detector and event selection cut efficiencies [11]. 
As stated above, at the lowest WIMP mass of 
6 GeV at which we show an exclusion limit in [2] , 
our result would only improve by a factor of 1.3, the 
improvement becoming less and less significant for 
higher WIMP masses. The reason why the limits 
improve is this: as the expected WIMP distribu- 
tion moves downwards in log 10 (S2/Sl), its separa- 
tion from the electronic recoil background distribu- 
tion increases, which in turn enhances the signal- 
to-background discrimination. The effect is fully 
exploited in the profile likelihood analysis, which 
does not use a pre-defined discrimination cut in the 
S2/S1 versus SI space. Nonetheless, such an analy- 
sis would not increase the sensitivity of the experi- 
ment significantly, while introducing additional sys- 
tematic errors. We have thus decided to publish a 
conservative result using a better known, Sl-based 
energy scale and to define the WIMP signal region 
based on 241 AmBe neutron calibration data down 
to the lowest WIMP masses that can be probed by 
our experiment. 

We add four additional comments to highlight the 
problems with the analysis performed by P. Sorensen: 

• The overall cut acceptance, as shown in Figure 
1 in [2], is not constant at 85%, as assumed in 
[1], but varies with energy, decreasing to ^50% at 
Sl=3 photoelectrons. In addition, the acceptance 
of the S2>150 photoelectrons condition is set to 
zero below Sl=l photoelectron, as shown in the 
same figure and detailed in [12]. This is effectively 
equivalent to not considering events below 3 keV 
nuclear recoil energy, where no £ ff measurements 
exist. Ignoring these cuts, as apparently done by P. 
Sorensen, will lead to an overestimated sensitivity 
increase, in particular at low WIMP masses. 

• The "detector-specific details" used in [1] were ob- 
tained from [13], which describes our detector dur- 
ing its commissioning run in 2009. The science 
data presented in [2] were acquired under different 
conditions, which will affect the analysis. For in- 
stance, the electron lifetime, determined by the liq- 
uid xenon purity, was considerably higher during 
the last run, an effect which has a direct impact 
on the size of the detected S2 signals. Moreover, 
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we use only the S2 signals detected by the bottom 
PMT array for computing the S2/Sl-ratio [2], and 
the electron extraction field was slightly lower in 
this run. 

• For our result reported in [2], the uncertainty in 
£ e ff is taken into account in the profile likelihood 
approach and a limit with 90% C.L. coverage is 
derived. This has to be included in the analysis 
presented in [1], otherwise the result is too opti- 
mistic. In addition, it is not clear if and how the 
known background expectation and its uncertainty 
for the 225 live-days run is taken into account. 

• We agree on the importance of understanding the 
population with low S2/Sl-values below the nuclear 
recoil region. A detailed study of these events, 
which are also present in our gamma calibration 
data of the electronic recoil band, is in progress. We 
nonetheless remark that the two observed events in 



XENONlOO's 225 live-days run are located within 
the bulk of the nuclear recoil distribution, and 
could thus be nuclear recoil candidates. At the 
same time, unlike the statement in [1], these events 
are also consistent with "anomalous leakage", as 
defined in [2] . 

To summarize, while we agree with P. Sorensen's con- 
siderations on the expected signal region for low-mass 
WIMPs [1] , we have observed only a mild impact on our 
analysis results, in contrast to the claimed "one order 
of magnitude or more" improvement in sensitivity. At 
present, we refrain from using this approach, for a robust 
analysis should be based on measurements rather than 
indirectly inferred quantities. Once direct measurements 
of the charge yield of nuclear recoils at low energies will 
become available, we will undoubtedly make use of this 
additional information to redefine our energy scale and 
strengthen our analysis potential for low-mass WIMPs. 
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